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Abstract 

We study systems of few two-component fermions interacting in a Harmonic Oscillator trap. 
The fermion-fermion interaction is generated in a finite basis with a unitary transformation of 
the exact two-body spectrum given by the Busch formula. The few-body Schrodinger equation is 
solved with the formalism of the No-Core Shell Model. We present results for a system of three 
fermions interacting at unitarity as well as for finite values of the S-wave scattering length 02 and 
effective range r2- We also present results for the four- and five-fermion systems at unitarity. In all 
cases, we observe a fast convergence of the many-body energies with respect to the size of the basis 
in which the effective two-fermion interaction is defined. The many-body energies are in excellent 
agreement with exact solutions for the three-body problem, and results obtained by other methods 
in the other cases. 
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I. INTRODUCTION 



Tremendous experimental progress has been made in the area of ultracold atomic sys- 
tems over the past decades. The abihty to dial the interactions among the particles via 
Feshbach resonances as well as the high-level of control over the external confinement, have 
enabled to utilize these systems as simulation tools for various phenomena, from condensed- 
matter physics to quantum information [l| . A recent ground-breaking achievement j2|, Isj is 
the ability to further confine just a few atoms in nearly isolated sites of optical lattices in 
which atoms are cooled down to extremely low temperatures below the millikelvin. In the 
limit of low tunneling, each lattice site may be regarded as an independent harmonic oscil- 
lator (HO) well. Theoretical models of many-body quantum systems that were considered 
rough idealizations in the past can now be realized experimentally and studied. By pro- 
viding well-defined, strongly interacting many-body systems, the study of ultracold atoms 
gives physicists an unprecedented opportunity to assess theoretical techniques that cross the 
boundaries of disciplines. 

In this paper, we explore one of these strongly interacting systems, namely an ensemble of 
two-component fermions confined in a HO potential. The temperature and density of matter 
reached in these ultracold dilute systems are such that the collisions among particles occur 
at small energy and the average interparticle distance is much larger than the range of the 
interaction, which in the case of neutral atoms, is given by the van der Waals length. As a 
consequence of this separation of scales, the interaction among the particles can be modeled 
by a zero-range potential, such as the pseudopotential 0] which depends on a few low-energy 
parameters, e.g., the S-wave two-body scattering length a2 and the effective range r2. For 
the two-particle system in a HO trap, analytical solutions with the pseudopotential exist 
and are given by the so-called Busch formula j5|-0]- Predictions of this exact model have 
been tested against experimental data and found to provide an accurate description of the 
two-atom system js], Isf. For systems with a larger number of particles, besides the three- 
body system at unitarity where analytical solutions exist joj, one has to rely on numerical 
calculations. 

In this work, we report on a new numerical study of few-fermion systems in a HO poten- 
tial. The solution of the few-body problem is obtained by solving the Schrodinger equation 
in a truncated HO basis using the formalism of the No-Core Shell Model (NCSM) [lo| . Since 
we work in a finite basis, an effective interaction has to be used instead of the pseudopotential 
which is only meaningful in the infinite Hilbert space. This effective two-body interaction is 
generated with a unitary transformation of the exact two-body spectrum given by the Busch 
formula. By construction the exact solutions of the two-particle system are reproduced in 
the truncated space. Unitary transformations have been widely used in different fields of 
physics. In Nuclear Physics effective nucleon-interactions have often been constructed using 



the Lee-Suzuki method [10Hl2|. Starting from a "bare" nucleon-nucleon potential defined 
in the infinite Hilbert space, the two-nucleon Schrodinger equation is first solved within 
a as-large-as-possible (but finite) two-body model space assumed to ressemble the infinite 
Hilbert space. The Lee-Suzuki transformation is then applied so that the interaction in the 
truncated space reproduces eigenenergies (usually the lowest) obtained with the bare inter- 
action in the large two-body model space. In this study of HO-trapped fermions, the unitary 
transformation will be directly applied to the exact analytical spectrum given by the Busch 
formula. Using this interaction, we will show that the energies of A-fermion systems (we 
will consider A=3,4,5) are in excellent agreement with results obtained by other methods. 
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Moreover, we will show that the convergence with respect to the size of the basis in which 
the effective interaction is constructed, is fast. 

Several other numerical studies of few-fermion systems in a HO trap have been previously 
reported in the literature. Different approaches based on Monte Carlo technique such as, 
the Green's function Monte Carlo [131, the fixed-node diffusion Monte Carlo UM and the 



Shell-Model Monte Carlo [15[ methods have been applied to these systems. In |16l-ll9|. the 



NCSM formalism has been used to study few-fermion systems in a HO trap with an effective 
interaction constructed as a series of contact interactions and derivatives of delta functions. 



according to the principles of Effective Field Theory |20H22| . Using the ground and few 
excited-states of the exact two-particle spectrum, a two-body effective interaction tailored 
to the truncation of the model space, is constructed order by order in a systematic manner. 



In [23[, the many-body problem is solved with the Configuration Interaction (CI) to study 
systems in the unitary limit with an effective interaction whose parameters are fixed such 
that the two-body eigenenergies given by this interaction correspond to the energies in the 
two-body exact spectrum. Our construction of the interaction has some similarity with that 
approach in the sense that, in both cases, all the energies of the two-body spectrum obtained 
with the effective interaction, correspond to the exact solutions given by the Busch formula. 
The authors in 23| assume a separable form for the two-body interaction whereas in our 
case, no such assumption is made. Moreover, we will consider not only the unitary limit but 
also finite values of the scattering length and the effective range. 

The paper is organized as follows. In Sec. Ullwe recall some generalities about the two- 
particle system in a HO trap interacting with a short-range interaction and we show the 
construction of the effective two-body interaction. In the following section. Sec. IIIH we 
show results for systems of three, four and five fermions and compare with results obtained 
by other methods. We conclude and discuss future applications in Sec. IIVI 



II. FORMALISM 



We consider a system of A identical fermions of mass m interacting in a HO trap of 
frequency uj. For definiteness, we think of two-component fermions but the framework can 
be straightforwardly applied to other fermions as well as bosons. The Hamiltonian Hiab for 
this system in the laboratory frame reads: 



Hlab = W— 



2 2 

muj r- 



i<j 



(1) 



where (p^) is the position (momentum) of particle i with respect to the origin of the HO 
potential (fixed in the laboratory) and Vij the two-body interaction between particles i and 
j. A nice property of the HO potential is the decomposition into a piece acting on the center 
of mass (CM) of the A particles and a piece acting on their relative coordinates. Removing 
the CM part, we obtain the Hamiltonian H describing the intrinsic motion of the particles 



H = Hn 



with 



1 

A 



E 



{Pi -fjf muj^ 
2m + 2 



(2) 



(3) 
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Using the formalism of the NCSM lOj, we will diagonalize the Hamiltonian if ([2]) in a finite 



basis constructed with the eigenfunctions of the HO potential. The two-body (effective) 
interaction Vij in ([2]) will be tailored to the truncation of the model space and constructed 
from the exact solutions of the two-particle system. 



A. Two-fermion effective interaction 

We begin first by recalling some generalities about systems interacting at low energy, and 
then we describe the construction of the fermion-fermion interaction for systems in a HO 
trap. 

It is well known that at sufficiently low energy, a physical system can be described without 
the complete knowledge of the interactions between its constituents. In free space, for 
a system of two particle with a relative momentum k and interacting via a short-range 
potential, the details of the interaction at short-distance r < R, with r being the interparticle 
separation and R the range, are irrelevant to describe the physics at /c <C 1/-R. Moreover, 
due to the presence of the centrifugal barrier, only the lowest partial waves contribute and 
the S-wave dominates. In this low-energy regime, the S-wave phase shift 6o{k) is given by 
the Effective Range Expansion (ERE): 

kcot6o{k) = -— + '^k^ + ..., (4) 

where a2, r2, ... are, respectively, the scattering length, effective range, and higher ERE 
parameters not shown explicitly. For energies close to the scattering threshold, the S-wave 
scattering length primarily determines the scattering of two particles, as it can be seen 
in Eq. (j4]) when /;;—)■ 0. Similarly, for a few(many)-body system at sufficiently low-energy, 
in which interparticle distances are large compared to the range of the interaction, the 
scattering properties will also depend mainly on jl^]. Physical properties of systems that 
depend only on the scattering length are universal. Universality in physics usually refers to 
situations in which systems that are very different at short distances have identical long- 
distance behavior. Two universal systems may have completely different internal structures 
and interactions, but the low-energy behavior will be the same as long as they have the same 
scattering length. Corrections coming from details of the interaction, such as the effective 
range r2 in Eq. (jlj) are called non universal. 

Ignoring short- distance details, it is then possible (and often practical) to substitute the 
original potential acting in the inner region, by the zero-range pseudopotential [3], 



^(5) _ 27r^^ tan(^o(A;)) d 



where /i is the reduced mass of the two-particle system. By construction, solutions of the 
Schrodinger equation with the pseudopotential (|5]) have the same asymptotic behavior as 
the original potential. The pseudopotential being of zero-range, the asymptotic behavior is 
reached for r — )■ whereas for the original potential, it is reached beyond distance of order 
of the range R. Hence, the wavefunctions given by both potentials differ only in the inner 
region. Equivalently, the solution of the Schrodinger equation with the pseudopotential ([5]) 
can be obtained by solving the free Schrodinger equation supplemented by the following 
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boundary condition at the origin, 



with C, a normahzation constant. 

We now consider a HO potential of frequency u in which two particles are trapped and 
interact via a short-range interaction. As in free space, we assume that the interaction 
between them can be replaced by the pseudopotential ([5]), or equivalently, by the boundary 
condition ([6]). The HO potential in the relative frame is given by V^^{r) = ^nuj^r"^. For 
r > 0, the S-wave wavefunction </)(r) of the two-particle system is solution of the Schrodinger 
equation 

\2i^d^ + 2'"^^'^^ '^^'^^ = ^'^^^^ > ^- 

At the origin, the HO potential V^^{r) vanishes and (j){r) fulfills the boundary condition 
(E]). Solutions of Eq. (JTj) can be written with the Tricomi confluent hypergeometric function 



U 25 



as 

„2 



0(r) = Are-^U \^^\^^ 3/2, rVfe^j (8) 

where A is a normalization constant and h = \ — the HO length. From the behavior of the 
Tricomi function near the origin, one obtains 

cx l^^^^^^lj^r - for r ^ 0. (9) 

where we use the notation v = By matching the behavior of 0(r) at the origin (see 

Eq. (Q) with the boundary condition ([6]), one obtains the Busch formula js]- 



3/A-E/2 \ 



/4-g/2 ^ 2 



kcotSoik). (10) 



The above equation relates the energy of the two particles in the HO trap to the S-wave 
phase shift 6o{k) in free space. Using the ERE (jlj), one can rewrite Eq. ffTOl) as: 



3/4-g/2 \ 

- ^ + (11) 



V4-i?/2 \ 2a2 



The Busch formula (fTT]) was initially demonstrated in (5| for interactions characterized with 
the ERE truncated to the scattering length 02 and was latter generalized to include more 
terms from the ERE j^, 0] . 

At this point, a remark should be made about the validity, for a trapped system, of 
replacing the original short-range potential by the pseudopotential (jS]) (or equivalently of 
using the boundary condition (^). As we wrote previously, the pseudopotential substitutes a 
short-range interaction by a zero-range interaction that gives the same asymptotic behavior. 
If the HO trap has a too large contribution in the inner region of the original potential, the 
matching with asymptotic continuum wavefunction looses its meaning since particles cannot 
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be considered free anymore. As a consequence, the HO should be small in this region, that 
is, b should be large enough i.e. b/R ^ 1 so that the matching procedure is not spoiled 
by the HO potential. A study on the limitation of using the pseudopotential for a trapped 
system can be found, for instance, in [6|. 

We now describe the construction of the effective two-body interaction. Using the eigenen- 
ergy solutions of the Busch formula fllOl) and the corresponding wavefunctions ([H]), one can 
write the Hamiltonian H^"^^ of the two-particle system in the trap as 

ff(2) ^ jj\E^^^U, (12) 

where E^'^^ is the diagonal matrix formed with the eigenenergies and U is the matrix formed 
with the associated eigenvectors^. The Hamiltonian H^'^^ is defined in the infinite Hilbert 
space. An obvious choice of a basis for this space is the set of relative S-wave HO states 
|n) characterized by the HO radial quantum number n and the energy En = {2n + 3/2)hu. 
Introducing the parameter nmax, we perform the following partition of the two-particle space: 
the HO states with radial quantum number smaller or equal to nmax define the model space 
V, within which the effective interaction will be constructed. The remaining states, i.e., 
the states with radial quantum larger than nmax, define the complementary space Q. The 
largest number of quanta N^'"^^, carried by states of V is N^^^ = 2nmax- Introducing the 
projectors P and Q acting respectively in V and Q, the matrix U in Eq. f|T2|) can be divided 
into four blocks: 



Upp UpQ 

Uqp Uqq 



(13) 



Each block is obtained by application of the projectors on U: for instance, by definition 
UpQ = PUQ. Similarly, E^^) can be split as 



where Epp {Eqq) is the diagonal matrix formed by the eigenenergies Eq. fllUp restricted to 
V {Q). The effective Hamiltonian H'^p is constructed such that the energies in the truncated 
model space, i.e. V, correspond to the lowest energies of the original Hamiltonian H^'^\ We 



construct with the same unitary transformation as in 26 1 



Hf = , "^^^ E?l , . (15) 

When nmax oo, the model space V converges to the full Hilbert space and the matrix 
product UppUpp becomes equal to the Identity operator. In this limit, the effective Hamil- 
tonian Hf^ converges to the original Hamiltonian H^^\ The two-body effective interaction 
Vp^ acting in P is obtained by subtracting the HO potential to the effective Hamiltonian 
Hf: 

Vf = Hf - PHoP. (16) 



For each eigenencrgy, the eigenvector can easily be expressed in any basis from the corresponding radial 
wavefunetion ([S]). 
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Thus, we have designed a two-body effective interaction Vp that reproduces the lowest 
energies of the two-particle system in the trap given by the Busch formula. From now on, we 
write Kifax interaction in the two-body space characterized by the cutoff A'^^'^^ = 2nniax- 

We briefly describe in the following, the construction of the NCSM basis for the resolution 
of the many-body Schrodinger equation ([2]). 



B. Many-body basis 



The resolution of the Schrodinger equation for the few-fermion system in a HO trap comes 
down now to the diagonalization of the effective Hamiltonian H^"^^^: 



i<j 



(17) 



where V^^^^ is the effective two-body interaction between particles i and j. The above 
Hamiltonian is diagonalized using the formalism of the NCSM where the basis is constructed 
from HO wavefunctions. Here we work with Jacobi coordinates [l^] defined in terms of 
differences between the CM positions of sub-clusters within the A-hodj system, e.g.. 




I A 



A-l 



A 



A 



(fi + f2 H h ta-i) - ta 



(18) 



The basis states are eigenstates of the HO potential Hq (j3]) expressed in Jacobi coordinates. 
For the purpose of illustration, let us consider a system of A=3 fermions of spin s = 1/2. 
In that case two Jacobi coordinates ^1,^2 are necessary to describe the intrinsic motion and 
the antisymmetrized basis states can be written as: 



A[[(t>nl{il) ® te(6)]^ \{SS)SS- S)} 



(19) 



where ra, / (A^, £) are respectively the radial quantum number and angular momentum as- 
sociated with the Jacobi coordinate ^1 {^2)- The angular momenta / and C are coupled to 
the total angular momentum L whereas the three spins s are coupled to total spin S. The 
antisymmetr y is enforced by the operator more details on the antisymmetrization can 
be found in 27|. The state f lT9|) is an eigenstate of the HO Hamiltonian with the energy 
(A^3 -|- 3)^0; where is the quantum number defined by iVa = 2?7, -|- / + 2M + L. For a 
system of A particles, the model space is truncated by introducing a cutoff A^™*^^ defined as 
the largest number of quanta in the eigenstates of HO used to construct the A-body basis. 
Using again the three-body system as a concrete example, truncating the three-body basis 
at N™^"^ means keeping only states with HO energies such that Nj, = 2n + l + 2N+L < N^^^. 

Instead of working with Jacobi coordinates, we could have chosen to construct the NCSM 
basis as a set of Slater determinants from single-particle HO wavefunctions defined in the 
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laboratory frame. The antisymmetry is easily built into this approach but the dimension of 
the basis space grows quickly with N^^^. On the contrary, working with Jacobi coordinates 
is more effective than a Slater-determinant basis in the sense that the dimension of the model 
space grows more smoothly with N^^^. On the other hand, with this choice of coordinates, 
the antisymmetrization becomes increasingly difficult as the number of particles grows. But, 
since in the current paper, we investigate systems with at most five fermions, we have chosen 
the Jacobi coordinates approach in order to be able to consider large A^^^^. 

We show in the following results for the energy of few-fermion systems in a HO trap for 
various values of the scattering length 02 and effective range r2. 



III. RESULTS 

We now present numerical results for systems of few two-component fermions interacting 
in a HO trap. The effective two-fermion interaction is constructed as described in the 
previous section. 



A. Three fermions at unitarity 

We begin here with an ensemble of three fermions at unitarity i.e. ^ = 00 and y = 
0. This systems gives us a perfect ground for testing our numerical approach since semi- 
analytical solutions exist in that case jof. In this regime, the only length scale characteristic 
of the system is given by the HO radius b. The energies of the two-fermion system at 
unitarity are solutions of the Busch formula ( fTTl) and are given by 

E = (1/2 + 2n)hu; (20) 

where n is integer such that n > 0. The two-body interaction VJ:^ is constructed in 
the truncated two-body model space characterized by the cutoff N^^^ = 2r2max and by 
construction, V^^^^ reproduces the nmax + 1 lowest energies in the spectrum fl20|) . We study 
the convergence of the three-body energy as the cutoff is increased. The three-body 

NCSM basis is characterized by the parameter N^^^ defined as the largest number of HO 
quanta carried by the states forming this basis. For each values of A^™^-^, the size of the 



three-body model space is increased until a convergence is reached [18|, [191, [23 

The ground-state solution of the three-fermion system is coupled to L'^ = 1~ and we 
show in Fig. [T]its energy for N^^^ = 10, 16 as a function of N^^^. For a fixed cutoff N^^"^, 
the Hamiltonian is fixed and one expects a variational behavior of the energy as the three- 
body model space is increased. This can be seen in Fig. [1] where for each value of N^'^^, 
the energy decreases sharply first as A'^^^^ increases, until a convergence to a final value is 
reached. The value of N^^^ for which the convergence is obtained depends on the particular 
value of iV™^^. For iV™*" = 10, the energy of the three-body ground state does not change by 
more than 10""^ when iV™^^ > 21 while for A'^^'^^ = 16 convergence at this level is achieved 
for N^^^ > 33 (see Fig. [T]). Figure [2] shows the convergence of the ground-state energy as a 
function of N^^^. The agreement with the exact solution is very good: for the largest cutoff 
we have considered here, A"™^^ = 20, we obtain 2.7712 for the ground-state energy (in units 
of hu) as compared to the exact solution 2.7727 [9|. Moreover, the convergence with respect 
to N^^^ is rather fast. Indeed, starting at N^^^ = 4, all the circles representing the energy 
of the ground state lie on top of the exact result represented by the dashed-line (see Fig. 
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max 

FIG. 1: Energy in units of the HO frequency, E/fiuj, of the ground-state L'^ = 1^ of the A = 3 
system at unitarity, as function of the three-body model-space size, N^^^: iV™'^-^ = 10 (square) 
and iV™'^^ = 16 (triangle). The dashed hne marks the exact value 



|2]). At 7V™^^ = 4, the energy obtained is 2.7657 which corresponds to a difference with the 
exact result smaller than 0.3 %. 

The calculation of excited-state energies can be carried out in a similar fashion. Again, 
for each value of N^^^, we increase N^^^ until convergence at 10~^ is reached. For the 
first-excited state coupled to L'^ = 1~, convergence at a fixed A^™*™ is reached for larger 
values of N!j^^^ as compared to the ground state: for N^^^ = 10 and N^'^^ = 16, one has to 
go to respectively to > 23 and A^3 > 35 to reach the same precision. Results for the first 
excited-state energy as a function of A^^^'*^ are shown in Fig. |2j As for the ground state, the 
agreement with the exact value is excellent: for N^^^ = 20, the energy is equal to 4.7732 
as compared to the exact value 4.7727 jof. Relatively to the ground state, the first-excited 
energy converges more slowly with N^'^^ (see Fig. [2]). But still, the energy approaches 
rapidly the exact solution: already at N^'"^^ = 6, the energy is 4.7769 which corresponds to 
a difference with the exact value 4.7727 smaller than 0.1 %. 

We have also studied states with other couplings. The general features of convergence 
with N^^^ and A^*^^ are similar as previously. We show in Fig. [31 the results for the ground 
state and first-excited state coupled to L'" = 0"*" as a function of A^™*". The convergence 
with respect to A^^^^ is fast and in excellent agreement with the exact values: for the ground 
(first-excited) state, at N^^^ = 22, the energy is equal to 3.1640 (5.1622) whereas the exact 
energy is 3.1662 (5.1662) @. 



B. Finite scattering length and effective range 

The construction of the two-body effective interaction is general and allows us to consider 
also physics away from unitarity. As one moves away from a Feshbach resonance, the effects 
of the interaction range should become more pronounced and universality reduced to some 
degree. Here, we will assume the following values for the scattering length and the effective 
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FIG. 2: Energy in units of the HO frequency, of the ground state (circle) and first excited 
state (diamond) coupled to L'^ = 1~ of the A = 3 system at unitarity, as function of the two-body 
cutoff jiV™'^^. The dashed lines mark the exact values for the ground state and the first excited 
state 
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FIG. 3: Energy in units of the HO frequency, E/Tiuj, of the ground state (circle)and first excited 
state (diamond) coupled to L"" = 0"*^ of the A = 3 system at unitarity, as function of the two-body 
cutoff, N^^^. The dashed lines mark the exact value 



range: b/ao = l,r2/b = 0.1. This choice is motivated by the fact that this system has been 

studied in [l8| using interactions constructed from the EFT principles, and allows us to test 

our method away from unitarity. 

Figure |4] shows the energy of the ground state and the first excited state coupled to 
= 1~ as a function of N^^^. As previously in the unitary regime, for each value of N^^^, 

the three-body model space is increased until convergence at 10~^ is reached. Again, the 
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convergence with N^^^ is fast: for the largest cutoff we have considered, N^^^ = 20, the 
energy for the ground state and the first excited state are respectively 1.8458 and 3.8739. 
Our results agree nicely with the energies obtained in the EFT approach: at N^LO, the 
ground-state and first-excited state energy were respectively found at 1.8443 and 3.8680 
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FIG. 4: Energy in units of the HO frequency, E/hu, of the ground state (circle)and first excited 
state (diamond) coupled to L'^ = 1~ of the A = 3 system for b/a2 = l,r2/h = 0.1 as function of 
the two-body cutoff, N2^''^^- The dotted line is drawn to guide the eye. 



C. Four- and five-fermion systems at unitarity 

We now consider heavier systems with A=4 and A=5 fermions. We show only results 
at unitarity although finite values for the scattering length and effective range could be 
considered as well. 

Again for each N^^^, the Schrodinger equation is solved in a truncated many-body model 
space whose size N^^^ is increased until convergence is reached. Because of the larger number 
of particles, we limit ourselves to smaller space. 

Figure [5] shows the ground-state and first-excited state energy of the four-body system 
coupled to L'^ = 0^ as a function of N^^^. The value A'^^"'^^ has been increased for each 
jymax yj^j^^[\ energies do not change by more than 10~^. For the largest cutoff considered 
here, N^^^ = 10, we obtain a ground-state energy equal to 3.537. This result is in good 
agreement with calculations using EFT interaction at N^LO where the energy was found to 



be 3.52 [18|, and also with results from various different approaches, where the ground-state 
energy was found to be 3.6 ± 0.1 3.551 ±0.009 and 3.545 ± 0.003 jlsj. For the 



first-excited state we obtain 5.085 which is also in good agreement with EFT results at N^LO 
where the energy was found to be 5.07 



Finally, we present results for the ground state energy of the five-fermion system at 
unitarity. Figure E] shows the ground-state energy for several cutoff N^^^, as a function of 
jymax^ j^^ prcvious cases, the energy at fixed N^^^ decreases when the many-body model 
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FIG. 5: Energy in units of the HO frequency, E/fioj, of the ground state (circle)and first excited 
state (diamond) coupled to L'^ = 0"^ of the A = 4 system at unitarity as function of the two-body 
cutoff, iVa"^''- 




2 4 6 8 10 12 14 



max 

N5 

FIG. 6: Energy in units of the HO frequency, E/hiv, of the ground state of the A = 5 system 
at unitarity, as function of the five-body model-space size, N^''^^: N2^^^ = 2 (square), N^^'^^ = 4 
(triangle) and N^"^"^ = 4 (star). 

space grows, and convergence is faster for smaller N^^^. For the largest five-fermion model 
space we consider here, i.e. for N^^^ = 13, results have converged nicely at N^^^ = 2 
and N^^^ = 4. For this latter cutoff, the difference between the energy at N^^^ = 11 and 
jymax — jg ^ 12 10~^, which represents less than 0.2 % of the energy 6.036 obtained at 
jymax _ p^j, jymax _ -j-j^^ difference in the energy for the same values of A'^^^^ is, as 
expected, larger. It is equal in that case to ~ 1.6 10~^ and this difference corresponds to 
2.6 % of the energy 6.068 obtained at N^^^ = 13. For all N^^^ we have considered, results 
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at N^^^ = 13 are all within a small range of energy that goes from 5.976 to 6.068 (5.976 
is the ground-state energy for N^^^ = 2). For the largest cutoffs N^^"^ = 4 and N^^^ = 6, 
the results for the ground-state energy are respectively 6.036 and 6.068. These values are in 
perfect agreement with the energy 6.1(1) obtained with the GFMC approach fl3| . 



IV. CONCLUSIONS AND OUTLOOK 

We have studied systems of few fermions in a HO trap interacting via a two-body effec- 
tive interaction generated from a unitary transformation of the exact two-body spectrum. 
Constructions of interactions tailored to a truncated space using unitary transformations 
are common in many fields of physics. The originality of our approach comes from the fact 
that the unitary transformation is directly performed on the exact two-body spectrum which 
is independent of the form of a initial "true" two-body potential. Indeed, the energies of 
the two-body system in the HO potential are given by the Busch formula as function of 
few model- independent, low-energy parameters, namely the S-wave scattering length a2 and 
effective range r2. This model independence is enabled by the separation of scales inherent 
to these systems. Indeed, the interparticle distance being much larger than the range of the 
interaction and the matter density much lower than the typical matter densities, the physics 
is only sensitive to the long-distance (low-energy) components of the interaction which can 
be taken into account, in a model- independent fashion, by the ERE. By construction, the 
two-body effective interaction reproduces, in the truncated space, the energies solution of 
the Busch formula. 

Using this effective two-body interaction, we have considered several few-fermion systems 
and studied the convergence of their energies with N^'^^ the interaction cutoff. The few- 
body Schrodinger equation has been solved using the formalism of the NCSM. We first 
studied systems of ^4=3 particles at unitarity and at finite values of 02, We have showed 
an excellent agreement with the exact solutions at unitarity for the ground state and the 
first-excited state. For finite values of 02 and we have seen a similar agreement with 



the results obtained by the EFT approach of [18|. In all cases, we have observed a fast 
convergence of the energies with respect to the cutoff N^^"^. We then considered systems of 
^4=4, 5 fermions at unitarity and similarly, we have obtained an excellen t ag reement for the 



energies with results obtained from several other numerical approaches [13|, [Ij, [181, [23 . 

The construction of the effective interaction from the Busch formula is quite general and 
can be used to study other systems such as few-fermion mixtures consisting of different 
species, and bosonic systems. An interesting application could be the study of the evolution 
of an ensemble of few bosons in a one- dimensional (ID) trap, from the weakly-interacting 
regime of Bose-Einstein Condensate (BEC) to the strongly-repulsive regime of the Tonks- 
Girardeau (TG) gas j35|. In this latter case, the infinite repulsion between bosons imposes an 
effective Pauli-principle, and the many-body wavefunction can be mapped to a many-body 



non-interacting fermion wavefunction. Several approaches |36H38| have been considered 
to study these ID systems. In the most recent one [11], the authors have developed an 
elegant method where the many-body wavefunction is constructed from a correlated pair 
wave function (CPWF). The CPWF is inspired from the idea that the analytical solution 
for a single pair in the HO trap can be extended to the many-body system by taking 
products of exact two-body functions. By construction, this method reproduces the exact 
known solutions for many-body systems in the limits of zero and infinite interactions. For 
the two-particle system, the exact value solution of the Busch formula for any arbitrary 



13 



interaction is reproduced. The CPWF method has been apphed to ID bosonic systems and 
two-component few-fermions mixtures and has shown an overall good agreement with other 



approaches, especially in regime close to the TG gas or weakly-interacting systems [38 . 

Our unitary-transformation-based approach has some similarities with the CPWF. In- 
deed, by construction, the two-body spectrum given by the Busch formula is also exactly 
reproduced in our case. It would then be interesting to compare results obtained by both 
methods for the momentum distributions, the density of the bosons and the occupation num- 
ber in the HO trap shells with the evolution of systems from the weakly-bound limit (BEG) 
to infinitely repulsive interaction (TG gas). Work in that direction is currently underway 
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